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Abstract. We investigate the kinetics of phase separation for a mixture of rod- 
like viruses (fd) and polymer (dextran), which effectively constitutes a system of 
attractive rods. This dispersion is quenched from a flow-induced fully nematic 
state into the region where the nematic and the isotropic phase coexist. We 
show experimental evidence that the kinetic pathway depends on the overall 
concentration. When the quench is made at high concentrations, the system 
is meta-stable and we observe typical nucleation-and-growth. For quenches at 
low concentration the system is unstable and the system undergoes a spinodal 
decomposition. At intermediate concentrations we see the transition between 
both demixing processes, where we locate the spinodal point. 



PACS numbers: 82.60.Nh,82.70.Dd,64.70.Md 

1. Introduction 

Systems that are quenched into a state where at least one order parameter is 
unstable undergo spinodal phase separation. Here, the initially homogeneous system is 
unstable against fluctuations of arbitrary small amplitude, and phase separation sets in 
immediately after a quench. In the initial stage of phase separation an interconnected 
" labyrinth structure" of regions with somewhat higher and lower values of the order 
parameter is observed. For quenches where the system becomes meta-stable, phase 
separation is initiated by fluctuations with a sufficiently large amplitude. Since such 
fluctuations have a small probability to occur, phase separation sets in after a certain 
delay time, referred to as the induction time. Here, nuclei are formed throughout the 
volume which grow when they are sufflcicntly large. The two different mechanisms 
of phase separation (spinodal decomposition and nucleation-and-growth) can thus be 
distinguished during the initial stages of phase separation from (i) the difference in 
morphology (interconnected structures versus growth of isolated nuclei) and (ii) the 
delay time before phase separation sets in (no delay time for spinodal decomposition 
and a finite induction time for nucleation-and-growth). As Onsager showed in 19490, 
the situation is different when the particles are not spherical in shape, i.e. disk-like or 
elongated particles. Here the system can become unstable or meta-stable with respect 
to fluctuations in orientation. These orientational fluctuations drive concentrations 
differences, resulting in a phase with high concentration and orientational order, the 
nematic phase, and a phase with low concentration and no orientational order, the 
isotropic phase. For very long and thin rods with short-ranged repulsive interactions, 
the binodal concentrations, i.e. the concentrations of the isotropic and nematic 
phases in equilibrium after phase separation is completed, have been determined using 
different approximations in minimizing Onsager's functional for the free energy (see 
ref. 12] and references therein) , while for shorter rods computer simulations have been 
performed to obtain binodal concentrations |31 0] . Also the spinodal concentration 
where the isotropic phase becomes unstable has been foundPE]. 

Binodal points are relatively easy to determine experimentally, since they are 
given by the concentrations of the bottom and top phase after phase separation. In 
contrast, it is not at all straightforward to obtain spinodal points, since one would 
ideally n like to perform a concentration quench from low or high concentration into 
the two-phase region, where the initial state is isotropic or nematic, respectively. In 
a recent paper such a kind of 'quench' was performed by inducing polymerization of 
short actin chains 6 , and factoids and spinodal structures were observed. Signatures 
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Figure 1. The bifurcation diagram, where the orientational order parameter 
P2 is plotted against concentration. Indicated are the various meta- or unstable 
regions for the two different initial states of the homogeneous suspension. The 
points marked by X and O are spinodal and binodal points, respectively. 



of spinodal decomposition have also been obtained for boehmite rods, by homogenizing 
a phase separated system and sequential polarization microscopy and Small Angle 
Light Scattering measurements||2|- Alternatively, external fields like shear flow 8, and 
a magnetic field |21 llOj can be applied to prevent a system from phase separation 
and to stabilize the nematic phase. After cessation of such an external field the 
nematic phase will become unstable or meta-stable, depending on the constitution of 
the sample, and phase separation sets in. In this paper we induce a fully nematic 
phase with a well defined director by imposing shear flow to a dispersion of colloidal 
rods. We use /rf-viruses as system, since the equilibrium phase behavior concerning 
the binodal points, has been well understood on the basis of Onsager theory [TTl I12j . 
Polymer is added to the dispersion in order to widen the region of isotropic-nematic 
phase coexistence, which facilitates the phase separation experiments 13 . We 
perform quenches of a flow aligned initial state to zero shear, which renders the 
system unstable or meta-stable to fluctuations in the orientation, depending on the 
concentration of rods. As a consequence phase separation sets in, which we observe 
by Confocal Laser Scanning Microscopy (CSLM). We perform this experiment for 
different concentrations, throughout the region of phase coexistence. Our results 
illustrate the difference between nuclcation-and-growth and spinodal decomposition 
in the case of demixing elongated particles, and result in the determination of the 
nematic- isotropic spinodal point. 



2. On the Instability of Initial States 

A convenient way to analyze the stability of a homogeneous initial state is to derive 
an equation of motion for the order parameter tensor, 

So(i) ^ j) duuuPoiu,t) . (1) 

The largest eigenvalue P2 of the tensor Qo{t) = | [ So(t) — 1 1] (where I is the identity) 
measures the degree of alignment. For the isotropic state P2 = 0, while for a perfectly 
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aligned state P2 — 1- A stability analysis of stationary solutions of this equation of 
motion is most conveniently made on the basis of a bifurcation diagram |S] , where P2 
for stationary solutions is plotted against the concentration. A schematic bifurcation 
diagram is given in fig^ The two solid lines represent stable stationary solutions of the 
equation of motion, while the dotted lines represent unstable stationary solutions. The 
isotropic state ceases to be stable above the concentration indicated as C^p'^"' , while 
the nematic state becomes unstable at concentrations lower than C^^^" . Above (7*^™ , 
the isotropic state is still a stationary solution , but is now unstable. Below C^''™, 
on the contrary, there is no unstable nematic state that is a stationary solution of the 
equation of motion. The two spinodal concentrations C^^^" and C*^™ are connected 
by a separatrix which separates the basins of attraction for the isotropic and nematic 
state. A homogeneous initial state above this separatrix develops a higher degree of 
alignment, while an initial state below the separatrix becomes more isotropic. 

Note that the bifurcation diagram relates to homogeneous systems. In 
an experiment, starting from a homogeneous state, inhomogeneities develop 
simultaneously with the change of the order parameter of the otherwise homogeneous 
system. In equilibrium, after completion of phase separation, there is an isotropic 
phase with concentration Cf™ in coexistence with a nematic phase with concentration 
C^'". One can either start from a stationary state, in which case Pq in eq.QJ is 
independent of time, or from a non-stationary state, like a nematic state with a 
concentration lower than C„, in which case the time dependence of the temporal 
evolution of alignment of the otherwise homogeneous system couples to the evolution 
of inhomogeneities through the time dependence of Pq . 

In this paper we prepare an initial nematic state shearing a suspension at large 
enough shear rate such that the induced nematic phase is stable against phase 
separation (see refs.|14| for a discussion of the bifurcation diagram for sheared 
systems), and then quench to zero shear rate. For this initial state it is expected that 
spinodal decomposition occurs at lower concentrations, while nucleation and growth is 
observed at higher concentrations. For an isotropic initial state this would be reversed 
: spinodal decomposition at high concentrations and nucleation and growth at lower 
concentrations. The observed phase separation kinetics thus depends crucially on the 
preparation of the initial state of the suspension. 

3. Materials and methods 

As model colloidal rods we use /d-virus particles which were grown as described 
elsewhere jl3_. A homogeneous solution of 22.0 mg/mL /d-virus and 12.1 mg/mL 
of Dextran (507 kd, Sigma- Aldrich) in 20toM tris buffer at pH 8.15 with lOOmM 
NaCl is allowed to macroscopically phase separate. This concentration of /d-virus is 
exactly in the biphasic region, which is very small when no polymer is added, namely 
between 21 and 23 mg/niL. Due to the added polymer, the binodal points shift to 
17 and 30 mg/mL, respectively. New dispersions are prepared by mixing a known 
volume of the coexisting isotropic and nematic bulk phases. The relative volume of 
nematic phase in this new dispersion is denoted as (pnem- 

For the microscopic observations we used a home-built counter rotating cone-plate 
shear cell, placed on top of a Leica TCS-SP2 inverted confocal microscope. This cell 
has a plane of zero velocity in which objects remain stationary with respect to the 
microscope while shearing. For details of the setup we refer to reference[^. For the 
measurements described here we used confocal reflection mode at a wavelength of 488 
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Table 1. Overview of the used samples. 



Code 


J' 


4 

^nern 






^nem 


fd {mg/niL) 


29.5 


28.1 


25.8 


23.6 


19.3 


^nern 


0.96 


0.85 


0.68 


0.52 


0.18 




Figure 2. The initial stages of phase separation for five different concentrations 
after a quench from a flow aligned nematic phase to zero shear. The top row 
shows the micrographs taken by reflection confocal scanning laser microscopy 
(field of view = 110 lim); the middle row shows the fourier transform of the 
micrographs; the bottom row plots the mean intensity of the micrographs minus 
the mean intensity for the isotropic phase, normalized by the initial intensity of 
the nematic phase. 



nm. Quench experiments were done as follows. Samples were first sheared at a high 
rate of 10 s^^ for several minutes. The shear was then suddenly stopped, after which 
images were recorded at regular time intervals. These images were parallel to the 
flow-vorticity plane. The table gives an overview of the concentrations where quench 
experiments have been performed. 



4. Results 



In the top row of Figure|21we show micrographs of the initial stage of phase separation 
for five different concentrations taken after a shear rate quench from a high shear rate, 
where the nematic state is stable for each concentration, to zero shear. These images 
show the flow (vertical) - vorticity (horizontal) plane at a given time after the quench. 
Thus the director of the initial nematic phase lies in the vertical direction. Fourier 
transforms of the images are plotted in the second row of Fig. |21 The background is 
corrected for by subtracting the Fourier transform of the first frame. The third row 
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Figure 3. The cross section of the Fourier transform parallel (a) and 
perpendicular (b) to the director for 'Pn^m- The length is scaled by the rod 
length L. 



plots the development of the total intensity of the images minus the intensity in the 
isotropic phase, as determined from an isolated isotropic region, normalized by the 
initial nematic intensity. Qualitatively the difference between the concentrations is 
obvious. In the first two images, i.e. the two highest concentrations, isolated dark 
ellipsoidal structures can be seen on a bright back ground. These are droplets of 
the isotropic phase referred to as factoids. The number of factoids increases when the 
concentration is decreased (b and c) until the structures become interconnected (d and 
e). This also follows from the fourier transform of the pictures where a ring is detected 
for the lowest concentration and a constant increasing intensity towards K=0 for the 
highest concentration. The time scale at which the inhomogeneities are formed also 
changes. As can be seen in the third row of Fig. |21 the high concentrations all show 
an induction time before the phase separation sets in, while for the low concentrations 
phase separation sets in immediately. Note also the times at which the images in Fig. 
[3 were taken. The isolated nuclei and the induction time are typical for nucleation- 
and-growth, while the interconnected structures and immediate phase separation are 
typical for spinodal decomposition. 

We use the Fourier transform of the images as shown in Fiq. |21to quantify the 
phase separation processes. The interesting quantity for nucleation-and-growth is the 
width of the Fourier transform, Afc, which is a measure for the anisotropic form factor 
of the nuclei. Alternatively one could determine the average size of the features in 
real space, but due to the low contrast this is difficult. For spinodal decomposition 
the interesting quantity is the wavevector at which the fourier transform reaches its 
maximum, k„iaxj quantifying the fastest growing concentration fluctuation. In both 
cases the fit of the Fourier transform should be performed in two dimensions, since the 
initial state is anisotropic. Therefore we took cross sections parallel and perpendicular 
to the director in the Fourier domain, i.e. the vertical and horizontal in Fig. |2] middle 
row, to determine kmax- Typical cross sections for are shown in Fig. O where the 
wavevector k is scaled by the rod length L. To determine A/c, we performed a 2-D 
gaussian fit around the origin of the Fourier transforms. Results of a 2-D gaussian fit 
of the Fourier transform around the origin are shown for the higher concentrations in 
Fig. 2Ji and b, plotting the width in the direction of the director and perpendicular 
to the director, respectively, kmax as found from fits of the cross sections parallel 
and perpendicular to the director are given in Fig. 01 c and d, respectively. Both fit 
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Figure 4. Width of the 2-D gaussian fit, AkL , of the Fourier transform parallel 
(a) and perpendicular (b) to the director for the higher concentrations for the 
higher concentrations. Wavevector kmaxL where the intensity is maximum for 
the cross sections parallel (c) and perpendicular (b) to the director for the lower 
concentrations. 
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Figure 5. (a) The ratio Afcxi/AfcyL for the higher concentrations and (b) the 
ratio fcmaa!,±^/'Smai, II ^ the lower concentrations. 



procedures result in an anisotropic morphology as can be seen in Fig. [SJ where we 
plotted Akj^L/AkiiL and kjnax,±L/k^ax,\\L- 

The late stages of the different phase separation processes also show some 
interesting phenomenology, as can be seen in Fig. El For spinodal decomposition we 
observe first a growing of the interconnected structures, which then break down into 
factoids. Later on factoids coalesce, and they become more spherical with increasing 
size. Note that these factoids contain the nematic phase and not the isotropic phase, 
as observed for the nucleation-and- growth process at higher concentrations. In the late 
stage of nucleation-and- growth, i.e. at high concentrations, we see that coalescence 
of factoids containing the isotropic phase as shown in the bottom row of Fig. is 
favorable when two factoids meet somewhat from the middle. In this case the rod 
orientation near both features is similar and the barrier which has to be overcome for 
coalescence is low. 
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Late stages spinodal decomposition 




coarsening spinodal structure splitting up 
t = 450s t= 1300 s t = 3800 s 




tactoid formation coalescence to spherical droplets 



Late stages nuoleation-and-growth 
f=500s t = 531s t= 563 s t = 594 s t= 625 s 




Figure 6. The late stages for spinodal decomposition in the top two rows (i^nemi 
field of view=375 /irri), and coalescence of tactoids in the bottom row (iPnemi held 
of view=73 fim). 



5. Discussion 

On the basis of these observations we can now locate the metastable region, i.e. where 
the system has to overcome a free energy barrier, and the unstable region, where there 
is no such barrier. At the high concentrations (<^^emi Vnem) the system is meta-stable, 
which is reflected by the observed isolated structures formed (top row in Fig[21) and 
the induction time (bottom row in Figl^J. With decreasing concentration the system 
approaches the unstable region: the number of nuclei increases while the induction 
time decreases and finally vanishes. The lowest concentration f^^jy^ is clearly unstable 
after cessation of the flow. It shows all the features typical for spinodal decomposition: 
phase separation immediately sets in throughout the whole sample, with a typical 
length scale which is characterized by the scattering ring observed in the Fourier 
transform. It can be shown, in fact, that the observed phase separation process for 
the lowest concentration has features typical for the spinodal decomposition of rods, 
as derived recently from a microscopic theory by one of the authors [T^. This will be 
the subject of a following paper [T7j. 

In the intermediate region it is difficult to judge from the morphology if 
nucleation-and-growth takes place or spinodal decomposition, since it is difficult 
to distinguish between a high number of tactoids and an interconnected structure. 
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However, <^^e^ shows a short induction time after the quench after which clearly 
separated tactoids are formed, while for (y5^g„ phase separation immediately sets in 
showing ellipsoidal structures which clearly 'influence' each other. Moreover, Figures 
0] shows that the size of the structures formed in (Pn^,^ coincides after some time with 
the clearly nucleated structures of (p^^^ and fnerm while the size of the structures 
formed in (Pnem coincides with samples which clearly show spinodal decomposition. 
Thus, we locate the transition from meta-stable to unstable, i.e. the spinodal point, 
between at 23.5 mg/mL and 25.8 mg/mL. This is the first experimental observation 
of the spinodal point in a rod-like system. We should mention at this point that in 
fact our sample consists of a mixture of rods and polymer. Addition of the polymer 
causes a widening of the biphasic region i-G- a shift of the binodal points. It 
is now interesting to see that the high concentration binodal shifts as much as from 
23 mg/mL to 30 mg/mL. In contrast, the high concentration binodal point, C^™, 
shifts from a concentration between 21 mg/mL and 23 mg/mL to somewhere between 
23.5 mg/mL and 25.8 mg/mL. This leads to the interesting conclusion that the shift 
of the high concentration binodal point, C^™, due to the attraction between the rods, 
is considerable compared to the shift of the high concentration spinodal point, C^^™. 
In other words, making the rods attractive causes a widening of the meta-stable region, 
while the unstable region remains unaffected. Addition of more polymer will result in 
more complex kinetics as described in reference |18|. 

Interestingly, for all concentrations we observe that the morphology of the phase 
separating system is anisotropic. This is most clear for the highest concentrations, 
where the tactoids all point upwards, i.e. in the direction of the director of the 
surrounding nematic phase. Also the Fourier transforms for the lower concentrations 
show deformed intensity rings in fourier space (most right FFT image in Fig. [^J. 
Moreover, also the kinetics of the phase separation is fastest in the direction of the 
nematic director. This follows for instance from the ratio of k^ax as plotted in[S|3, 
which increases in time. In other words, for all concentrations phase separation is 
anisotropic, due to residual alignment after the quench of the initially strongly sheared 
suspension, and not isotropic as is the case for spheres [T^. 

The length of the first observed tactoids just below C^™ is about 12 times the 
rod length, while just above C^^™ it is seven times the rod length. The thickness 
is about two third of the length in both cases. Typical length scales for the initial 
spinodal morphology are not more than six rod lengths. These sizes seem to be quite 
small, considering also the random orientation of the rods in the isotropic phase, 
but it is in accordance with a the microscopic theory for spinodal decomposition of 
rods^B]. It does suggest that we really image the initial stage. The breaking up of 
the spinodal structure into nematic tactoids and the sequential growth in the late 
stage of spinodal decomposition seems surprising since for dispersions of spheres only 
coalescence and macroscopic phase separation would be observed. However, a similar 
order of events has been observed for polymer mixtures with thermotropic liquid 
crystals 1201 • Simulations on such mixtures show that the break down is due to the 
effect of the flow- alignment coupling, and not primarily due to elastic effects j^J. An 
explanation in the same line was given by Fukuda in a numerical treatment of time- 
dependent Ginzburg-Landau equations of liquid crystalline polymers [221 ■ The volume 
dependence of the morphology in the final stage can be explained by the competition 
between the interfacial tension and nematic elasticity of the tactoids [2^1 ■ 
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6. Conclusion 

We studied the kinetics of the nematic-isotropic phase transition of a dispersion of fd- 
virus particles with added polymer after shear quenches into the two-phase region. By 
varying the equilibrium rod concentration ipnem we were able to detect a nucleation- 
and-growth mechanism for high ifnem, spinodal decomposition for low ifnerm and the 
transition between the two processes. In this way we were able to trace for the first 
time the nematic-isotropic spinodal point C^^™. Thus, we found that addition of 
polymer widens the meta-stable region greatly. Furthermore, we showed that the 
phase separation is strongly influenced by the director of the initial nematic state. 
The nematic phase also influences the late stages of spinodal decomposition, causing 
a splitting up of the interconnected structures. 
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